Quantum Capacities of Bosonic Channels 
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We investigate the capacity of bosonic quantum channels for the transmission of quantum in- 
formation. Achievable rates are determined from measurable moments of the channel by showing 
that every channel can asymptotically simulate a Gaussian channel which is characterized by second 
moments of the initial channel. We calculate the quantum capacity for a class of Gaussian channels, 
including channels describing optical fibers with photon losses, by proving that Gaussian encodings 
are optimal. Along the way we provide a complete characterization of degradable Gaussian channels 
and those arising from teleportation protocols. 



One of the aims of Quantum Information Theory [lj is 
to follow the ideas of Shannon and to establish a theory 
of information based on the rules of quantum mechanics. 
A key problem along this way is the calculation of the 
quantum capacity of noisy quantum channels. That is, 
the question how much quantum information — measured 
in number of qubits — can be transmitted reliably per use 
of a given channel? Despite substantial progress |2| this 
can be answered only in very few cases Q as a simple 
formula, comparable to Shannon's coding theorem, is not 
known. 

In this work we investigate the quantum capacity of 
bosonic channels which might describe transmission in 
space, as light sent through optical fibers, or in time, like 
in quantum memories |^|. The paper has two parts. In 
the first part we prove that the quantum capacity of any 
bosonic channel T is lower bounded by that of a cor- 
responding Gaussian channel Tq, which can be derived 
from measurable moments of T . This implies that for 
determining and certifying achievable rates for the trans- 
mission of quantum information through T we need not 
know the channel exactly (which might be hardly possible 
in infinite dimensions), but merely its second moments, 
i.e., a few measurable parameters. In the second part 
we then explicitly calculate the quantum capacity of a 
class of Gaussian channels, which includes the important 
case of attenuation channels modelling optical fibers with 
photon losses and broad-band channels where losses and 
photon number constraints might be frequency depen- 
dent. Along the way we provide two tools that might be 
of independent interest: a complete characterization of 
degradable Gaussian channels and of those arising from 
teleporting through Gaussian states. 



PRELIMINARIES 

Before we derive the main results we will briefly recall 
the basic notions 0, Q . Consider a bosonic system of N 
modes characterized by N pairs of canonical operators 
(Qi, Pi, ■ ■ ■ , Qn, Pn) ='■ R for which the commutation 



relations [Rk,Ri] — io~ki are governed by the symplec- 



p2N 



tic matrix a. The exponentials We 
are called Weyl displacement operators. Their expecta- 
tion value, the characteristic function, x(£) := tr [pl 7 ^] 
is the Fourier transform of the Wigner function and for 
Gaussian states it is a Gaussian 



X(0 = e 



(1) 



with first moments d k = tr [pRk] and covariance matrix 
(CM) T ki := tr [p{R k - d k ,R t - di}+\. Note that coher- 
ent, squeezed and thermal states in quantum optics are 
all Gaussian states. 

Gaussian channels 0, transform Weyl operators as 
h- > Wx£,e~~^ Y ^ and act on covariance matrices as 



7 1 



(2) 



Particularly important instances of single-mode Gaussian 
channels are attenuation and amplification channels for 
which X = y/fjl and Y = \r) - 1|1. For < 77 < 1 this 
models a single mode of an optical fiber with transmis- 
sivity 77 where the environment is assumed to be in the 
vacuum state. The latter reflects the fact that thermal 
photons with optical frequencies are negligible at room 
temperature. For 77 > 1 the channel becomes an ampli- 
fication channel, where the noise term Y is now a conse- 
quence of the Heisenberg uncertainty. 

Teleportation channels: We will now derive the 
form of Gaussian channels which are obtained when tele- 
porting through a centered bipartite Gaussian state. As 
this is useful for applying but not necessary for under- 
standing the following it might be skipped by the reader. 

Let r = ( ^ c ) be the CM of a Gaussian state of 

Na + Nb modes with Na = Nb- Assume Bob wants 
to teleport a quantum state of Nb modes with CM 7 to 
Alice. Using the standard protocol he sends pairs of 
modes from 7 and Tb through 50:50 beam-splitters, mea- 
sures the Q and P quadratures, and then communicates 
the outcomes. Depending on the latter Alice applies dis- 
placements to the modes in T^. The simplest way of 
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deriving an expression for the output is to start with the 
Wigner representation and to assume that the state to be 
teleported is a centered Gaussian. The Wigner function 
before the measurement is up to normalization given by 
exp — £,[Mg S (T 7)~ 1 A/bs]£, where Mbs corresponds 
to the beam-splitter operation. With £ = 
the final Wigner function is then proportional to 



-£ [mImI s (r© 7 ) - 1 MbsMx] t 



(3) 



where Mx incorporates the displacements, i.e., it is the 
identity matrix plus an arbitrary 2Nb X 2Nb off-diagonal 
block which maps the 2Nb measurement outcomes onto 
the respective displacements. In order to circumvent in- 
tegrating Eq. we can now go to the characteristic func- 
tion, i.e., the Fourier transformed picture. The integra- 
tion then boils down to picking out the upper left block 
of the inverted matrix [M^M^ S (T ® j)- 1 M B sM x ]' 1 ■ 
The inversion is, however, trivial since Mgg = Mg S and 
M^ 1 is obtained from Mx by changing the sign of all 
off-diagonal entries. In this way we obtain that the input 
CM is transformed to 



X T -yX 



r. 4 + T C AX + (T C AX) T (4) 
+X T A T T B AX 



where A = diag(l, — 1, 1, — 1, . . .) and X is such that \[2X 
is the matrix of displacement transformations, i.e., the 
gain which is typically chosen to be v21- 

Clearly, Eq.(0} has the form J2J an d following the 
above lines it is straight forward to show that the chan- 
nel is Gaussian and maps any (not necessarily centered 
Gaussian) input characteristic function %in into 

X«*(0 = Xin (XOxrtf e a*o . (5) 

For standard protocols (X = 1) on single modes (Na — 
Nb = 1) this was derived in ||. 



ACHIEVABLE RATES FOR ARBITRARY 
CHANNELS 

The subject of interest is the quantum capacity Q(T) 
of an arbitrary — a priori unknown — channel T. We will 
show how one can certify achievable rates for the trans- 
mission of quantum information through T by only look- 
ing at the CM T of a state pr = (T ® id)(V0 which is 
obtained by sending half of an arbitrary entangled state 
tp through the channel. T could be determined by ho- 
modyne measurements. The argument combines (i) the 
relation between entanglement distillation and quantum 
capacities observed in llfl . (ii) the extremality of Gaus- 
sian states shown in |ll| and (iii) the explicit form of 
Gaussian teleportation channels derived in the previous 
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Bob 



FIG. 1: In order to obtain a Gaussian channel from an ar- 
bitrary quantum channel T Bob (the sender) prepares n in- 
stances of an entangled state xp half of which he sends through 
T® n . After applying two arrays of 50:50 beam-splitters to the 
output p^ n = [(T <g> id)(^)]® n the n reduced states will con- 
verge to a Gaussian state Q(pr) (with the same CM as pr) 
which can in turn be used to establish a Gaussian teleporta- 
tion channel Ta- 



section. All together this leads to the chain of inequalities 



Q(T) > D^(p T ) > D^{Q{ PT )) > Q(T G ) 



(6) 



Here -D,_(/?t) is the distillable entanglement under pro- 
tocols with one-way communication (from Bob to Alice). 
Since a classical side channel does not increase Q(T) this 
is clearly a lower bound to the capacity as Alice and 
Bob could simply first distill pr and then use the ob- 
tained maximally entangled states for teleportation [Icj . 
The second inequality uses that replacing px by a Gaus- 
sian state G(pr) with the same CM V can only decrease 
the distillable entanglement ^lj. Finally, if we use the 
Gaussian state in turn as a resource for establishing a 
teleportation channel Tq we end up with the sought in- 
equality Q(T) > Q(Tq). Tq is then the Gaussian channel 
in Eas. (|4l5() . which is for a fixed teleportation protocol 
(a fixed matrix X) completely determined by T. 

Bounds on the quantum capacity of Gaussian channels 
were derived in |l2l Il3l | and we will show below that it 
can be calculated exactly for some important cases. Note 
that a simple bound for Q(T) can be obtained from a 
lower bound to D*-(Q(pt)), the conditional entropy of 
the Gaussian state with CM T, i.e., Q(T) > S(T A )~S{T). 

Before we proceed, two comments on the quality of 
the above bound and its operational meaning are in 
order: The given argument holds for arbitrary T and 
■ip. However, since we bound by Gaussian quantities 
the inequality might become trivial (i.e., Q(Tg) = 
though Q(T) ^> 0) if both T and ip are too far from 
being Gaussian. On the other hand, if T is Gaussian 
and \ip) — (coshr) -1 ^ n (tanhr)™|rm) is a two-mode 
squeezed state, then in the limit r — * oo the inequality 
becomes tight, i.e., Q{Tq) — * Q(T) with exponentially 
vanishing gap. 

The first and last step in Eq.© have a simple opera- 
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tional meaning: sender and receiver first establish some 
entanglement, distill it and then use it as a resource for 
teleportation. The second step can also be understood in 
operational terms. To achieve p? <— > G{pt) both sender 
and receiver have to apply an array of beam-splitters (see 
FigDJ to many copies of the shared state px- Asymptoti- 
cally every reduced state at the output will then converge 
to Q{pr) [ill ]. By applying this gaussification to many 
different subsets, sender and receiver can then distill from 
or teleport through independent copies of G(pt) and in 
this way asymptotically simulate the channel Tq via T. 



QUANTUM CAPACITY OF GAUSSIAN 
CHANNELS 

It was proven in [2j that the quantum capacity of a 
quantum channel T can be expressed as 

Q(T) = lim - sup j(p,T® n ) , (7) 

n— >oo n p 

J(p,T) = S(T(p))-S((T®id)ty)) , (8) 

where ip is a purification of p and J is known as the 
coherent information. In general, the calculation of Q(T) 
from the above formula is a daunting task since (i) the 
coherent information is known to be not additive, i.e., the 
regularization n — > oo is necessary, and (ii) due to lacking 
concavity properties there are local maxima which are 
not global ones. On top of this, for bosonic channels the 
optimization is over an an infinite dimensional space. 

Fortunately, for a class of Gaussian channels includ- 
ing the important case of the lossy channel, these ob- 
stacles can be circumvented by exploiting recent results 
on degradability of channels |3J, LLJ] and extremality of 
Gaussian states |Tl) . 

To this end consider a channel T(ps) — tr e[U(ps <8> 
<Pe)U^] expressed in terms of a unitary coupling be- 
tween the system S and the environment E which is 
initially in a pure state ifE- The conjugate channel 
T c (ps) = trs[U(ps®<pE)U'] is defined as a mapping from 
the system to the environment. As shown in 3] the coher- 
ent information can be expressed in terms of a conditional 
entropy if there exists a channel T" such that T' oT = T c 
in this case T is called degradable. More precisely, if 
pgi e' is the extension of the state ps> —T'o T(p) to the 
environment E' of T", then 

J(p, T) = S(ps'E') - S{ps>) =■ S(E'\S') . (9) 

The conditional entropy S(E'\S') is known to be 
strongly sub-additive i.e., for a composite system 
S(E' 12 \S' 12 ) < S^si) + S{E' 2 \S' 2 ). This has impor- 
tant consequences: for a set {Ti} of degradable channels 
J(p, ®iTi) < *}2iJ{pi,Ti), where pi are the correspond- 
ing reduced states, and if each Ti is a Gaussian channel, 



Q 




FIG. 2: Quantum capacity of a channel with photon losses 
as a function of the transmission length I in terms of the 
absorbtion length l a , i.e., n = e~ l ^ la . For quantum memories 
I and l a are storage and decay time. The capacity vanishes 
for I /l a = In 2 ss 0.693, where the channel can be considered 
to be part of a symmetric approximate cloning channel. 

we have in addition 

J{p,®iTi) < J(pi,Ti) <^J(G(pi),TJ. (10) 

i i 

The last inequality follows from the extremality of Gaus- 
sian states w.r.t. the conditional entropy together 
with the fact that for Gaussian channels T c can be chosen 
to be Gaussian and the CM is transformed irrespective 
of whether the input was Gaussian or not. As a con- 
sequence, if Ti are degradable Gaussian channels, then 

Q(®iTi) =J2 sup J(p G ,Ti) , (11) 

i PG 

where the suprcmum is now taken only over Gaussian in- 
put states pc- Calculating the latter for Gaussian chan- 
nels is now a feasible task which was solved for the single- 
mode case in and in for broadband channels 
under power constraints using Lagrange multipliers. In 
fact, if we impose a constraint on the input energy of 
the form J^. uJiNi = £, where is the average input 
photon number of mode i with corresponding frequency 
u>i, then the above argumentation still holds, since the 
constraint just depends on the CM. The importance of 
Ea. (|ll|) stems from the fact that a large class of Gaus- 
sian channels is indeed degradable, as shown in 0] and 
extended below. In particular we can apply Eci. (|ll|) to 
attenuation (amplification) channels with transmissivity 
•q (gain ,Jrf). Together with the optimization carried out 
in |l2( this yields Q as a function of 77 (see Fig. 2): 

Q(r ? ) = max{0 ) log 2 |r ? |-log 2 |l-j 7 |} . (12) 

Note that the quantum capacity of every degradable 
Gaussian channel can easily be calculated as J becomes 
a concave function of the CM such that local maxima are 
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global ones. This is again a direct consequence of Eq.JHl 
together with the concavity of the conditional entropy [lj 
and the extremality of Gaussian states 0, . 



DEGRADABLE GAUSSIAN CHANNELS 

We will now investigate the condition under which 
EcL ljllJI was derived and characterize the set of degrad- 
able Gaussian channels — extending the results of |LJ] . To 
this end we represent the channel in terms of a unitary 
coupling between the system with Ng modes and a (min- 
imally represented) environment of Ne < 2Ng modes 
which arc initially in the vacuum state with CM 7_e = 1. 
The interaction is described by a symplectic matrix of 
size 2(Ne + Ng) x 2(Ne + Ng) which we write in block 

form as S = y ^ ^ ^ 0] . The output CM of the channel 

T : 7 1— * X~{X T + Y is then simply the lower right block 
of S(j E © l)S T (i.e., D = X and Cj E C T = Y) whereas 
the conjugate channel T c corresponds to the upper left 
block. 

Let us first focus on the case Ng — Ne and assume 
for simplicity that the blocks in S are non-singular. A 
channel is degradable if T c o T -1 is completely positive 
which is for a Gaussian trace preserving map equivalent 
to the condition Q 



Y + icr > iXaX 1 



(13) 



Inserting the above block structure and using 15] shows 
that complete positivity of T c o T _1 is equivalent to 



< (1 + ia) - K{l + i<r)K T 
K = C T D- T (jD- x C . 



(14) 



Expressing this in terms of X and Y finally gives [lfil ] 

(2XaX T a T - l)Y > . (15) 

Similarly we can derive a condition for degradability of T c 
(anti- degradability of T) which is again given by the ex- 
pressions (|14l I15|) which have then to be negative instead 
of positive semi-definite. 

Since for Ne = Ng = 1 I is a 2 x 2 matrix and 
thus XaX T a T = JdetX, condition (|15(l implies that 
either T or T c is degradable, as shown in [14j . Hence, 
as anti-degradable channels have zero quantum capacity 
(due to the no-cloning theorem), the quantum capacity 
of every Gaussian channel with Ng — Ne — 1 can easily 
be calculated. In fact, by utilizing the freedom of act- 
ing unitarily before and after the channel (which does 
not change its capacity) one can bring the channel to a 
normal form [l7| which only depends on the symplectic 
invariant detX such that Q(T) of every such channel is 
given by Ea. (|12|l with r\ = detX. 

Let us finally briefly comment on the case Ne 7^ Ng. 
If the environment is smaller than the system, then we 



can easily follow the above lines for instance by choosing 
a re pres entation of the channel with larger Ne equal to 
Ng It is worth mentioning that if S corresponds 

to a passive (i.e., number preserving) operation, then for 
Ne < Ng there are always unaffected modes such that 
Q(T) — 00 without additional constraints. If Ne > Ng 
then Ea. (|15fl is merely a necessary, whereas Ea. (|14fl is 
still a necessary and sufficient condition for degradability 
[lij ] . Applying the latter to a general single- mode channel 
with Ng = 1,Ne — 2 shows that generically one has 
neither degradability nor anti-degradability. Hence, it 
remains open whether in this case the capacity is given by 
Ea. pi(l . However, we can easily derive an upper bound 
by exploiting the fact that every Gaussian channel T can 
be decomposed as T = T± o T2, where T2 is a minimal 
noise channel for which Ne — Ng with X2 = X, 
Y% < Y and T\ is a classical noise channel for which 
X\ = 1, Y\ = Y — Yo. Due to the bottleneck- inequality 
for capacities (cf. |l2|]) we have Q(T) < Q(T2) where the 
latter is in the single-mode case again given by Ea. H12(l 
with 77 = detX. A lower bound is always given by the 
r.h.s. of Ea. (jTT|) as calculated in [l^ . 

In summary we characterized the set of degradable 
Gaussian channels and showed that their quantum ca- 
pacity can be calculated as it is attained for Gaussian 
product inputs. For arbitrary non-Gaussian channels we 
derived a certifiable Gaussian lower bound. Both ideas 
can be applied to finite dimensional systems as well. This 
will, however, be the content of future work 
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